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1 Introduction 



Open systems typically give rise to resonances. A resonance is a long-living quasi- 
stationary state, which eventually decays into the continuum. Physically, it may 
be thought of as a particle, initially trapped inside the system, which eventually 
escapes to infinity. 

One common approach to studying resonances is based on the analytic properties 
of the scattering matrix S{£) in the complex energy plane. Resonances correspond 
to poles 

£n = — (1) 

of S{£) on the non-physical sheet [HE]. In an alternative equivalent approach, which 
we shall follow here, one solves the Schrodinger equation subjected to the boundary 
condition of purely outgoing wave outside the range of the potential. This boundary 
condition, which describes a process in which a particle is ejected from the system, 
renders the problem non-Hermitian. The Schrodinger equation with this boundary 
condition leads to complex eigenvalues which correspond to resonances [H [2] . For 
a recent lucid discussion of resonances in quantum systems, with particular emphasis 
on the latter approach, see [31 H]. 

The outgoing-wave approach leads, in a natural way, to non-Hermitian effective 
hamiltonians , whose complex eigenvalues are the resonances of the studied system 
[3 El [7]. Such effective hamiltonians are very useful for studying resonances in 
scattering theory, including scattering in chaotic and disordered systems [HI El [IHl 

There are many examples of resonances in atomic and nuclear physics. Recently, 
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there has been considerable interest in resonances which arise in chaotic and disor- 
dered systems. See [9] for a recent review. One of the main goals in these studies is 
computation of the distribution P{T) of resonance widths. There is ample amount of 
work on computing P(T) in one-dimensional disordered chains[T3l HH HH |16l ITTl [HI 
[T9] . Numerical results presented in some of these works indicate that P{T) ~ r~'^ 
in a large range of values of F, where the exponent 7 is very close to 1. 

A more general quantity than P{T) is th^ density of resonances (DOR) 

p{x,y) = '^5{x -Re£n)5{y -ImSn). (2) 

n 

It is widely believed that the averaged DOR in the complex plane contains infor- 
mation about the Anderson transition P, [TTl [121 [22] • This expectation is based 
on an analogy with Thouless' arguments concerning the sensitivity of eigenstates 
to the boundary conditions in Hermitian localization theory [20l |2T]. Indeed, the 
coupling of the disordered system to the external world plays in our case a role 
similar to changing the boundary conditions in Thouless' picture. Namely, the 
width of a typical resonance in the insulating regime should be exponentially small, 
^typ ~ exp —L/C,{E) {L being the size of the system), whereas in the metallic regime 
the typical width is Ttyp ~ V/L"^, namely, the inverse Thouless time scale (P is the 
diffusion coefficient in the disordered metal). Thus, Ttyp, measured in units of level 
spacing A, is analogous to the Thouless conductance. This picture was already 
pursued numerically in [22] . 

The continuum limit of the disordered chain was studied in [23]. For simplicity, a 

chain opened only at one end was studied. The spectral determinant for the problem 

^In order to avoid cluttering of our formulas, we do not use the resonance width r„ in as 
an argument of p, but rather y = — r„/2. 

learned this argument about the expected scaling behavior of the DOR from B. Shapiro. 
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was derived, and the averaged DOR was expressed in terms of a certain integral over 
the solution of a certain singular two-dimensional Fokker-Planck equation. (That 
Fokker-Planck equation determined the probability distribution of the logarithmic 
derivative of the outgoing wave at the open end of the chain.) 

The present work was motivated in part by [131 [H]. particular, an analytical 
approach was developed in [H] for studying resonances, which is based on counting 
poles of the resolvent of the non-Hermitian tight-binding effective hamiltonian of the 
open chain. In the case of a semi-infinite disordered chain, coupled to a semi-infinite 
perfect lead, these authors have derived an exact integral representation for the 
DOR, valid for arbitrary disorder and chain-lead coupling strength. In the limit of 
weak chain-lead coupling (in which resonances are typically narrow) they were able 
to rigorously derive a universal scaling formula for the DOR, valid for any degree 
of disorder and everywhere inside the unperturbed energy band of the closed chain. 
The 1 /r behavior of the DOR follows from that formula. 

In this paper we shall review and explain how to construct energy dependent 
non-hermitian hamiltonians for studying resonance statistics in open systems. While 
many (but by no means all) of the results presented in this paper are known, we be- 
lieve our presentation offers a somewhat fresh look at these issues. Upon elimination 
of the leads, one can reformulate the problem in terms of an effective non-hermitian 
hamiltonian, which depends only on the degrees of freedom of the disordered system. 
In this effective description, the outgoing-wave boundary condition in the original 
system is translated into a local non-hermitian, energy dependent boundary condi- 
tion at the contact points (or more generally, contact regions) of the system and the 
leads. 
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This paper is organized as follows. In Section 2 we discuss resonances in a 
generic quantum system coupled to the external world by a single one- dimensional 
lead ( a single channel lead). We derive a general expression for the DOR in terms 
of an appropriate diagonal matrix element of the resolvent of the original closed 
system. From this expression, we derive an integral representation for the averaged 
DOR of the disordered system. In Section 3 we specialize to the case of an open 
one dimensional disordered chain, derive the corresponding effective hamiltonian, 
and obtain its continuum limit. The resulting continuum effective non-hermitian 
hamiltonian differs from the hermitian one of the closed system by a complex energy 
dependent boundary condition. The structure revealed in this way is quite generic, 
and we conclude in Section 4 by mentioning similar continuum effective hamiltonians 
for higher dimensional systems. 

2 Resonances in a System Connected to a Single 
Perfect Lead 

In order to keep the discussion as simple as possible, let us consider a quantum 
system connected to a single perfect semi-infinite one-dimensional lead, which lies 
along the negative x-axis. This construction is described in Figure [Tj We shall 
model the lead by means of a tight-binding hopping hamiltonian, with nearest- 
neighbor hopping amplitude t. The sites on the lead lie at the points Xn = na, n = 
0, — 1, — 2, . . ., a being the lattice spacing. Let us assume that the (closed) quantum 
system has an A^-dimensional state space, and that it is described by an x 
hermitian matrix Hij, H = W . We further assume that the system lives on some 
graph with nodes, and that Hij is the matrix element connecting site i to site j 
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(the link < ij > is directed, of course). Let us now connect the lead's end n = to 
some site in the system, which with no loss of generality we pick to be site i = 1. 
The hopping amplitude along the buffer link < 01 > is taken to be t', which need 
not be equal to t. In particular, t' = corresponds to a closed hermitian system. 
Disorder is modeled by some probability distribution for the matrix H, which means, 
in general, both random hopping and random site energies on the graph. 




Figure 1: A disordered system is opened up and coupled to the external world 
by a perfect lead, stretched along the negative x-axis. The lead is modeled by a 
tight-binding hamiltonian, with nearest neighbor hopping amplitude t. The < 01 > 
link connects the lead and the system, with corresponding hopping amplitude t'. 
The Schrodinger equation is to be solved with outgoing wave boundary condition. 
In this figure we have set the lattice spacing a = 1. See the text for more details. 
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The Schrodinger equation for this system is therefore 





= zipn , n < -1 






= zlpo 5 n = 






= zipi n = 1 






= zipn , n > 2 . 


(3) 



Here we lumped the wave- function amphtudes inside the system into an A^- dimensional 
vector tp, and z is the complex eigenvalue. Imposing outgoing-wave boundary con- 
dition in the perfect lead means that 

^„ = Ae-*'="^ n<0, (4) 

where the wave-vector k must be restricted to the right half of the fundamental 
Brillouine zone, namely, 

< Re fca < TT , (5) 

as the wave propagates freely to the left, into the lead. (This choice has the obvious 
continuum limit Refc > describing free propagation to the left.) 

We shall now eliminate the lead entirely from (|3]), following the idea presented 
in [8]. To this end we first substitute Q in the first equation in (|3]), from which we 
find that 

z = z{k) = —2tcos{ka) . (6) 

Next, we eliminate ipo = A = je^'^'^ipi from the n = equation, and substitute it in 
the n = 1 equation. In this way we discover that 

- ^e^'^^V^i + (HiJ), = z{k)i^, . (7) 
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The remaining equations {Hip)„ = z{k)ijjn {n > 2), as well as the equation ([7|) for 
n = 1, can be neatly written as 

H - —e^''P\i) = z{k)il) , (8) 

where P = |1)(1| is the projector on site 1, to which the lead is connected. Thus, we 
can ignore the lead from now on, and describe the open system itself by an effective 
hamiltonian 

H,ff = H-7]e'''-P, (9) 

with 

V = j- (10) 

Hgff is a non-hermitian N x N matrix, which also depends on energy, through k, 
according to This explicit k dependence is of course, the price we had to pay 
in order to eliminate the lead from the description of our system, and it should 
actually be expected of an effective description - the effective hamiltonian normally 
depends on the energy scale one studies. Of course, when t' = rj = 0, the lead is 
disconnected, and H^jf = H of the closed system. 

The picture we have in mind is that starting with the closed system (77 = 0), all 
N energy eigenstates are real and sharply defined. Then we open the system adia- 
batically, i.e., increase 1] slowly. As a result, each sharply defined energy eigenstate 
in the original system should broaden continuously into a resonance, with complex 
energy z{k,r]). Thus, we end up with N complex resonance energies, the solutions 
of Hf,ff{k)ilj = —2t cos{ka) ip . 

Let us briefly elaborate on the domain in the complex-/c plane which corresponds 
to resonances. We shall follow the discussion in [21 S], which give probabilistic 



interpretation to the modulus square of the time-dependent resonance eigenfunction 
For this, we shall temporarily include the lead in our discussion. In a 
resonance (i.e., quasi- stationary) state the amplitude \l/n(t) = ?/'„e~*^*^'^^* has to grow 
in magnitude into the lead, where the particle is likely to be after a long time. For the 
same reason, it also has to decay as function of time, at any finite fixed site. Thus, 
in addition to (|5]), we must also demand that both Imk and Imz = 2tsm(Re ka) ■ 
sinh(Im ka) be negative. Thus, sin(Re ka) > 0, which holds automatically due to 
As was originally discussed in [5] (and later extended in we see that as 
time goes by, we can maintain the numerical value of the spatial integral of |\I^„(t)p 
(the probability) , provided we allow the integration domain to expand at a constant 
speed to the left, which is nothing but the ballistic velocity of the ejected particle. 
To summarize, resonances must all lie in the strip 

0<Refca<7r, lmfc<0 (11) 

in the fourth quadrant of the complex k plane. Similarly, anti-resonances, which 
describe a situation in which the system absorbs particles from the leacQ must all 
lie in the strip — tt < Re fca < , Im /c < in the third quadrant of the complex k 
plane. 

2.1 The Secular Equation 

Resonances are the roots of the equation 

det{z-Heff) = 0, (12) 

•^Here the wave should propagate freely in the lead to the right, towards the system. Hence 
—IT < Heka < 0. Furthermore, Imfc < 0, since at t = the it is overwhelmingly probable to find 
the particle in the lead, while at the same time we must also have Imz > 0, since the probability 
to find the particle at n = must grow. 
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with k lying in the appropriate strip (11). Practically, it is easier to compute the 
ratio of determinants 



where 



is the resolvent of H. In (12)-(14) we must of course set z according to ([6]). 

Note that {GP)nm = Gni^mi- Hence, 1 + rje^^"'GP is a lower diagonal matrix. 



and computation of the last determinant in (13) is immediate. We find simply that 



F{k) = l + rie""'Gii{z{k)) . (15) 

Thus, in order to solve for the resonance spectrum of our model, all we require is 
the Gil element of the Green's function of the original closed system. The latter is 
the Green's function of a hermitian hamiltonian, and therefore well-studied. Note 
that we have not specified the specific nature of the closed system corresponding to 
H . Our discussion is completely generic! 

2.2 The DOR 

For a given realization of if, F{k) is a holomorphic function of k, and has zeros at 
the eigenvalues of H^fj and poles (on the real axis) at the eigenvalues of H. Let fc^ 
and A;^ be, respectively, the zeros and (purely real) poles of F{k). Thus, 



^ h-h^ ^ h-hZ' ^^^^ 



From the identity 



F'{k) 

m ^k-ki ^k-ky 

d 1 



dk* k — q 
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7rS^^\k-q), (17) 



which is nothing but Gauss' Law in 2d electrostatics (for a unit point charge located 
at position k = q), we thus find 



Averaging this equation with its complex-conjugate, we finally obtain that 



(18) 



1 



log 



m =Y.^^'\k-ki)~Y.^^^\k-ki). (19) 



2-K dkdk* 

a p 

Since the poles live entirely on the real axis, going off it and into the fourth quadrant 
in the complex fc-plane, we obtain our desired DOR. 



Continuing the analogy with 2d electrostatics [21], observe that (19) is noth- 
ing but the Poisson equation, relating the charge distribution on the LHS, to the 
electrostatic potential 

2 

det(2; - Heff) 



W{k,k*) 



log 



F{k) 



log 



det(;z - H) 



(20) 



det(2 - H, 



on the RHS. Moreover, note that the real quantity 

portional to the determinant of the 2N x 2N hermitian operator 

z - Heff 



in (20) is pro- 



n 



Hlff 



(21) 



2;' - 

In fact, given a non-hermitian operator, such as H^ff, whose spectrum we wish to 
study, the method of hermitization [25] instructs us to construct its hermitized form 



(21), and study its spectrum, which of course lies entirely on the real axis. Thus, for 



example, the Green's function — Ti) is analytic in the complex ( plane, save for 
poles (or a cut, upon averaging) along the real axis, where the spectrum is located. 
Thus, one may bring the power of analytic function theory to bear in analyzing the 
spectrum, which cannot be done for the non-hermitian H^ff. 
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2.3 The Averaged DOR 



In our closed disordered system 

Gn{z) =X{z) + zY{z) (22) 
is a complex valued random variable, with probability distribution 

P(X, Y; z) = {5{X- X{z)) 5 {Y - Y{z))) , (23) 



which we assume to be known. Thus, from (19) and (23) we immediately obtain an 



integral representation for the averaged DOR as 



1 d' 



2 



2 



Pa.(fc,A;*) = ^^^ I dXdYV{X,Y-z{k))\ogl + r,e'^'^{X + iY) . (24) 



3 The One-Dimensional Disordered Chain and the 
Continuum Limit of its i/e// 

We shall now depart from the general discussion and take H to be the tight-binding 
hamiltonian of a disordered chain with N sites, i.e., the one-dimensional Anderson 
model. We take the nearest-neighbor hopping amplitudes to be t, as in the lead. 
The site energies e„ (n = 1,2,..., A^) are i.i.d. random variables taken from some 
probability distribution q{e). Thus, the corresponding hermitian matrix if in ([3]) 
(and in Figjl]), in the previous section, is given by 

Hmn = -t{6m,n+l + 5m+l,n) + Cn^mn , I <m,n < N . (25) 

The resulting Schrodinger equation is therefore 

- t{lpn+l + 'ipn-l) + enlpn = Zlpn , (26) 
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with Dirichlet boundary conditions 



i^o = i'N+i = , 



(27) 



corresponding to a closed chain. As can be seen from (|9j), the effective Schrodinger 
equation for the open system (with the lead ehminated, of course) is obtained from 



(26) (or (25)) simply by replacing e„ by 



e„ = e„ - rie''"'Sni ■ 



(28) 



Statistics of resonances in this model was studied in detail in [T3l [HI [15 



The effective Schrodinger equation {H^ff — z)ip = can be formally obtained by 
applying the variational principle 

5S 



5^^ 



(29) 



to the complex action 

N 



n=l 



a \ a a 



N 



-D aip. 



n=2 



" a? ' 
(30) 



where 



D = ta^ and (, = z + 2t = At sin^ 



ka 



(31) 



are, respectively, the diffusion constant (the lattice version of y~) ^"^^ the shifted 



(renormalized) energy, and 



(32) 



is the symmetric second difference (the lattice discretized version of the self-adjoint 
laplacian) . 
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Strictly speaking, we should really apply the variational principle to the real 



action corresponding to the hermitized form (21). However, in order to keep the 
discussion as brief as possible, and since all we want to obtain in this section is the 
continuum limit of H^ff, and not to pursue the averaged DOR in detail, we shall 
contend ourselves with the complex action S. 

3.1 The Continuum Limit 

The continuum limit is obtained by sending a — )■ and t — )■ oo simultaneously, while 
keeping D = to? = ^ and k finite. Furthermore, t' — oo as well, such that the 
ratio 

./\ 2 

-J =e'^ (33) 

with A finite. In this limit we also obtain the familiar relation ( = Dk^. In the limit, 
the lattice amplitudes tend to the continuous wave function, ipn = 'ipina) — )■ ipi^x) 
and the site energies tend to the potential e„ — )■ V{x) . Obviously, — )■ d1il){x) 



and ^^^^^ — )■ 'il)'{a+). Finally, of course, ^ — )■ 5{x — a) and Ylm^ ~^ I (2'''' 
tend to the continuum Dirichlet boundary conditions 

m = HL) = (34) 
with L = Na (where of course — t- oo). Plugging all these limiting quantities in 



(30), we obtain the continuum limit of S as 

L 

S = j dx [{V{x) - |^(x)|2 - Di)*{x)dli){x) - D5{x - a) [i)* {x)d^i){x) + (A + ik)\i){x)\^)] 

0+ 
L 

dxr{x){H!'}f-C)H^), (35) 

0+ 
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with the boundary conditions (34) understood. We can immediately read off the 



continuum effective effective hamiltonian from (35), namely, 

H^f}' = Dp'^ + V{x) - D6{x -a){zp + X + ik) (36) 
with p = —idx as usual. 



Note that we have explicitly left the infinitesimal lattice spacing a in (35) and 



(36) as a mnemonic. In fact, S{x — a) in these expressions really stands for a thin 
boundary layer around the left end of the chain, with a very large coefficient, which 
penalizes for having dxip{x) + (A + ik)il){x) 7^ in the immediate vicinity of x = 0+. 
The continuum Schrodinger equation derived by applying the variational principle 



to (35) generates in this way the continuum resonance boundary condition 

a,7/^(0+) + (A + zA;)7/;(0) =0, (37) 

which depends on energy, through k. Since in the lead, a; < we have the outgoing 
wave il){x) = ilj{0)e~'^^^ , the derivative dxip{x) jumps: iIj'{0+) — ipiO—) = — A^/'(0). 



This jump is the result of the singular contact potential term —XD5{x) in (36). Note 



from (36) that A —00 penalizes for having ip{0) 7^ 0. Thus, this limit corresponds 



to Dirichlet boundary conditions, namely, t' = and a closed chain, as can be seen 



also from (33) 



More precisely, for a very small but finite, one integrates the Schrodinger equation 
—Dd1il){x) + V{x)il){x) = C'^(x), subjected to i'iL) = 0, from the right end of the 
system all the way to x = a, where the large coefficient of the boundary layer 
interaction takes over, and fixes ^/''(a— ) = —(A + ik)il){a). The wave function has 
then to relax to zero at a; = across the thin boundary layer, with tremendous 
slope. This segment of the wave function is an artifact, which we cut and throw, 
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and replace by the resonance boundary condition (37). 



4 Concluding Remarks Concerning Higher Dimen- 
sional Systems 

The structure revealed by analyzing the continuum limit of the one dimensional case 
is quite generic. The effective non-hermitian hamiltonian is generically given by the 
original differential expression (in the coordinate representation), supplemented by 
an appropriate energy dependent complex boundary condition. In conclusion, let 
us mention briefly two simple 3d examples in the continuum, which correspond to 
having infinitely many weak channels connecting the system to the environment. 
These results are straightforward: 

• Disordered half-space coupled uniformly to the environment through a contact 
plane. 

Let's take the disordered system to live in the z > half-space, and let it 
communicate with its environment through the xy plane. For a given complex 
energy z = ^^^^ there is a continuum of resonances indexed by the components 
of the wave- vector k_|_, perpendicular to the z axis, which are real. They 
correspond to the direction in which the particle is ejected from the system. 
Let q = k^he the on-shell complex component of momentum in the z-direction, 
such that q"^ = — k'j_. Then, the resonance amplitude immediately outside 
the system, must satisfy the outgoing boundary condition 




(38) 
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We can then obtain dzipkj_{z = 0^) right inside the system by considering the 



contact potential —D\S{z), in complete analogy with (37). Since the boundary 



condition (38) is rotationally symmetric with respect to the 2;-axis, we expect 



the averaged DOR to inherit this symmetry as well. 

A disordered ball of radius a coupled uniformly to the environment through 
its surface. 

In this case, we should consider resonances with definite angular momentum 
quantum numbers /, m in the outside world. For a given complex energy 
z = ^^2^' corresponding outgoing wave amplitude ipim{r) = Aimhi{Qr) 
must be proportional to a Hankel function. Thus, it must trivially satisfy 

h'(u) 

^Lia+) = Qj^^iM , (39) 

hi{u) 

where u = Qa. Again, the radial derivative immediately inside the ball may 
be obtained by taking into account the jump in the radial derivative due 
to a uniform radial-shell contact potential. Due to spherical symmetry, the 



boundary condition (39) is independent of m. Consequently, the averaged 



DOR should inherit this property as well. 
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